Creep motion in a random-field Ising model 
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We analyze numerically a moving interface in the random-field Ising model which is driven by a 
magnetic field. Without thermal fluctuations the system displays a depinning phase transition, i.e., 
the interface is pinned below a certain critical value of the driving field. For finite temperatures 
the interface moves even for driving fields below the critical value. In this so-called creep regime 
the dependence of the interface velocity on the temperature is expected to obey an Arrhenius law. 
We investigate the details of this Arrhenius behavior in two and three dimensions and compare our 
results with predictions obtained from renormalization group approaches. 
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I. INTRODUCTION 



In recent years the understanding of driven interfaces 
has improved considerably. Well known models of such 
interfaces are the equations of Edwards and Wilkinson |l]] 
as well as of Kardar, Parisi, and Zhang g]. Of par- 
ticular interest are driven interfaces moving through a 
quenched disordered medium which exhibit a so-called 
depinning phase transition. Without disorder the veloc- 
ity of a driven interface grows linearly with the applied 
driving force or driving field, respectively. This behavior 
changes in the presence of quenched disorder. For small 
driving fields the interface is pinned by the disorder. The 
interface moves only if the driving field exceeds a criti- 
cal value, i.e., on increasing the driving field a continu- 
ous phase transition from a pinned to a moving interface 
takes place (see, for instance, || and references therein). 
The expected dependence of the interface velocity on the 
driving field is sketched in Fig. |l|. For very large driving 
fields the disorder can be neglected and consequently the 
velocity depends linearly on the driving field. The depin- 
ning transition happens due to the competition between 
the disorder and the driving field. The disorder induces 
some effective energy barriers that suppress the interface 
motion. The driving field reduces these energy barriers 
but they are overcome only if the driving field exceeds 
the critical value. Examples of systems exhibiting a de- 
pinning transition are charge density waves QJ^] , or field 
driven domain walls in ferromagnets || . 

However, in the above scenarios thermal fluctuations 
are neglected. In real systems these fluctuations occur 
and no critical behavior is observed for finite tempera- 
tures. The reason is that even below the critical driving 
field the energy barriers can be overcome due to thermal 
fluctuations, resulting in a moving interface. A strik- 
ing effect of thermal fluctuations occurs at the critical 
field where the interface velocity v depends on the tem- 
perature T according to v ~ T 1 ^ with an exponent 
5 > 1 Another effect of thermal fluctuations is 

the so-called creep motion that occurs for driving fields H 
well below the critical threshold at sufficiently low tem- 
peratures. Here, the interface velocity is expected to be 
characterized by an Arrhenius behavior 



H 



FIG. 1. Schematic sketch of the interface velocity v of the 
pinning phase transition and its dependence on the driving 
field H. The bold line corresponds to zero temperature, 
T = 0. For small but finite temperatures the critical behavior 
is smeared out (thin solid line). The creep regime for small 
driving fields is characterized by an Arrhenius like behavior 
with an effective energy barrier E. 



-E(H)/T 



(1) 



with a certain field dependent energy barrier E(H). 
Creep motion was investigated, for instance, within the 
theory of flux creep phenomena |lfj[ ] and in several renor- 
malization group approaches to the Edwards- Wilkinson 
equation fil| , |l2|| . Experimentally, the Arrhenius behav- 
ior of the creep regime was observed for magnetic do- 
main wall motion in thin films composed of Co and Pt 
layers @. 

In this paper we consider the interface motion occur- 
ring in a driven random-field Ising model (RFIM) in the 
creep regime. In the next section we describe the details 
of the model and the simulations. In Sees. [lll| and |y| we 
investigate numerically the creep motion of the interface 
in the two and three dimensional RFIM. In particular, 
we show that the velocity behavior can be described by 
an Arrhenius ansatz and we investigate the temperature 
and field dependence of the prefactor of the Arrhenius 
law and study the energy barrier E(H). In section [v] we 
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summarize and discuss our results. 



III. CREEP MOTION IN THE TWO 
DIMENSIONAL RFIM MODEL 



II. MODEL AND SIMULATIONS 

To study the creep regime we consider the RFIM on 
a square or simple cubic lattice of linear size L. The 
Hamiltonian of the RFIM is given by 



H 



H Si - h i Si 



(2) 



where the first term characterizes the exchange inter- 
action of neighboring spins (Si = ±1). The sum is 
taken over all pairs of nearest neighbor spins. The spins 
are coupled to a homogeneous driving field H and to 
quenched random-fields hi which we choose to be uncor- 
rected ((hihj) oc 5ij) with (hi) — 0. Throughout this pa- 
per we consider uniformly distributed disorder, i.e., the 
probability density p that the random field takes some 
value hi is given by 



p(hi) 



(2A)- 




for 



\hi\ < A 
otherwise. 



(3) 



Using antiperiodic boundary conditions an interface is 
induced in the system which can be driven by the field 
H (see H for details). Starting with an initially flat in- 
terface we apply a Glauber dynamics with random se- 
quential update and heat-bath transition probabilities 
(see, for instance, H]). In our simulations the interface 
moves along the [lipand [11 Indirection of a simple cu- 
bic lattice. This is a natural choice since in the absence 
of disorder interface motion occurs for any finite driv- 
ing field ||. This property is an advantage especially in 
the creep regime where the interface is driven by small 
driving fields at low temperatures. 

The basic quantity in our simulations is the velocity 
of the moving interface, which is determined in the fol- 
lowing way. The interface movement corresponds to an 
increasing magnetization which is monitored as a func- 
tion of time, i.e., the number of Monte Carlo steps per 
spin. Starting from a flat interface the system after a 
certain transient regime reaches a steady state where the 
average magnetization grows linearly in time. The veloc- 
ity of the interface is defined as the time derivative of this 
magnetization. Spin flips outside the interface may also 
occur, caused by the finite temperature. These isolated, 
rare spin flips are unstable for sufficiently small tempera- 
tures, i.e., they flip back in the next update. Thus these 
spin flips do not affect the measurement of the global av- 
erage magnetization time dependence and therefore do 
not affect the determination of the interface velocity. 

During its motion the originally flat interface roughens 
due to the disorder. The width of the interface increases 
and finally reaches a stationary state. For the data pre- 
sented in this paper we have verified that the interface 
width remains small as compared to the extension of the 
system perpendicular to the interface. 



We measured the velocity of the interface in the creep 
motion regime. Since the creep regime is "far away" from 
the critical point, we expect that finite-size effects can be 
neglected. Investigations of various system sizes confirm 
this assumption and wc use therefore in our simulations 
a large number of update steps instead of large system 
sizes. We performed for each temperature and field value 
at least 10 5 Monte Carlo steps. Additionally we focus our 
analysis on one value of the disorder strength (A = 1.2). 
We have also performed some spot checks at different 
values of A in order to confirm that the results are not 
sensitive to the disorder strength. 

As mentioned above the velocity is expected to obey 
an Arrhcnius law 



v(H,T) = C(H,T)e~ E W/ T 



(4) 



in the creep regime. The effective energy barrier E(H) 
is independent of the temperature and tends to zero for 
H — > H c . Following a renormalization group analysis pd| ] 
we assume that the temperature dependence of the pref- 
actor of the Arrhenius law is characterized by a power-law 
behavior 



C(H,T) = c(H)T~ 



(5) 



with some particular exponent x. Independent of the 
actual value of x the interface motion stops for any finite 
value of the energy barrier (H < H c ) in the limit T — > 0. 
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FIG. 2. The interface velocity v as a function of the 
temperature T for various values of the driving field 
(H e {0.25,0.3, 0.35, 0.6} from bottom to top) and 
A = 1.2. According to Eq. (|) we plot \nvT x vs 1/T. On 
varying the exponent x we obtained nearly straight lines for 
x — 0.89 ± 0.17. The cutoffs at low and high temperatures 
are caused by different effects. At low temperatures the in- 
terface can be pinned for finite time intervals depending on 
the particular disorder configuration. In this case the inter- 
face displays a stop and go behavior which results in strong 
velocity fluctuations (not shown for clarity). The cutoff for 
large v occurs because the creep regime is exited at high tem- 
peratures or high driving fields. 
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FIG. 3. The energy barrier E and the prefactor c [see 
Eqs. (^) and (^)] as a function of the driving field .ff. To 
avoid an overlap between the two curves we plot c — 1 instead 
of c. 



In the first step of our analysis we determine the ex- 
ponent x. In an Arrhenius plot hivT x vs 1/T the expo- 
nent x is varied until straight lines are obtained. Good 
results are found for x — 0.89 ± 0.17 and the correspond- 
ing curves are shown in Fig. ^. 

A regression analysis of these curves then yields the 
value of the prefactor c(H) and the value of the energy 
barrier E(H) [Eqs. (^) and (||)]. The results are plotted in 
Fig. |^. On increasing the driving field the effective energy 
barrier decreases as expected. But the prefactor of the 
Arrhenius law displays no significant field dependence, 
i.e., c(H) = const. This is confirmed by Fig. |^ which was 
obtained by plotting lnvT x as a function of E(H)/T. 
The curves for different values of the driving field are 
seen to coalesce to a single curve which happens only if 
c(H) = const. 

We analyze the field dependence of the energy barrier 
starting from a recently performed renormalization group 
approach assuming 



E{H) = E 



*V-i 

H 



(6) 



On approaching the critical field of the depinning tran- 
sition H c the energy barrier vanishes. The value of the 
critical field H c = 1.12 ± 0.02 is obtained from an inde- 
pendent simulation at zero temperature. Thus we plot- 
ted the rescaled velocities as a function of (H c /Hy L — 1 
and tried to obtain a coalescence of the data similar to 
Fig. ^ by varying the exponent \x. (See Figj|). Our anal- 
ysis shows that only a logarithmic field dependence for 
/i — > fits the data, i.e., 



E(H) = E„ In 



(7) 



(see Fig. ^[). As one can see in Fig. ^ the logarithmic 
Ansatz yields a quite convincing fit. Our analysis there- 
fore suggests that the effective energy barrier displays a 
logarithmic field dependence. 
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FIG. 4. The rescaled interface velocity \nvT x vs E(H)/T. 
The values of the energy barrier E{H) are obtained from a 
regression analysis of the corresponding curves of Fig. § The 
data collapse confirms that the prefactor c displays no signif- 
icant field dependence [Eqs. (^) and (|^)]. 



IV. THE THREE DIMENSIONAL MODEL 



We next analyze the interface velocity of the three 
dimensional model for A = 1.7. For this value of the 
disorder the critical behavior has been investigated and 
the corresponding critical field has been found to be 
H c = 1.37±0.01 M. As for the two dimensional model we 
have also performed some simulations for values different 
from A = 1.7 to ensure that the main results discussed 
below do not depend on the particular choice of A. 

By driving the interface at finite temperatures and 
fields below the critical threshold we measured the ve- 
locity v(H,T). Again, we fitted the simulation data to 
Eqs. (|i|) and (||) by varying x to get straight lines in 
the \nvT x vs 1/T plot. A good fit is obtained using 
x = 0.79 ±0.09 (Fig. @). 

The result of the regression analysis for E(H) and c(H) 
is shown in Fig. [§. As in the two dimensional case E(H) 
decreases with increasing driving field and the prefactor 
is essentially independent of the driving field. On plot- 
ting the rescaled interface velocities vs E{H)/T the data 
coalesce (see Fig. |9|) only if c(H) = const. 

We consider now the field dependence of the energy 
barrier. In analogy to the previous section we check the 
conjectured field dependence obtained from a renormal- 
ization group approach. Applying the data of the inter- 
face velocities to the Ansatz Eq. (p|) yields a similar result 
as in two dimensions, i.e., the accuracy of the data col- 
lapse increases for /i — > 0. Therefore we again assume 
that the dependence of the energy barrier on the driving 
field displays a logarithmic behavior [Eq. (Q)]. The cor- 
responding curves are shown in Fig. QOl As can be seen 
the logarithmic ansatz yields a good fit of the velocity 
data in the creep regime. 

But we have to admit that in contrast to the two di- 
mensional case we observe here that different Ansdtze for 
the field dependence of the energy barrier may also lead 
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FIG. 5. The rescaled interface velocity \nvT x as a function 
of [(H C /HY - 1]/T [see Eq. @]. No coalescence of the data 
could be obtained for any finite value of /x. But with decreas- 
ing exponent u the coalescence becomes better. The figure 
shows the corresponding curves for [i = 0.05. 
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FIG. 6. The rescaled interface velocity In vT x as a function 
of In (H c /H)/T. This logarithmic Ansatz yields a quite con- 
vincing data collapse, i.e., this result suggests that the field 
dependence of the energy barrier is given by Eq. (^) in the 
limit (i —¥ 0. 



to a fit of the data. For example, several authors conjec- 
tured that the energy barrier is given by 



E(H) = E 



He 
H 



(8) 



for H <^ H c . Note that Eq. (||) agrees with the one dis- 
cussed above [Eq. <Q)] for sufficiently small driving fields. 
The above Ansatz was derived within a theory of flux 
creep behavior jl0| and is expected to hold for the present 
situation of driven interfaces. In this case the exponent 
[i is given by fi — (2( + D — 3)/(2 — Q with £ denot- 
ing the roughness exponent of the interface at the depin- 
ning transition (see fl2|,[l3| and references therein). For 
the Edwards- Wilkinson equation with quenched disorder, 
£ = (5 — D)/Z has been determined by an e expansion 
within a renormalization group scheme fl5f| . This value 
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7. Interface velocities of 

the three dimensional model for different temperatures and 
driving fields (H — 0.3, 0.35, 0.4, 0.6, from bottom to top). 
On varying the exponent x we obtained nearly straight lines 
for x = 0.79 ± 0.09. 
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FIG. 8. The energy barrier E and the prefactor c [see 
Eq. (JE])] as a function of the driving field H. 



is believed to be exact to all orders of e [[L6| , and inserting 
it into the formula above yields /i = 1. independent of D. 

Fittin g ou r data according to Eq. (g) yields fi = 0.825± 
0.1 (Fig. [Fy). The accuracy of the fit is similar to the one 
obtained from the logarithmic Ansatz of E{H). Thus, in 
the three dimensional case one cannot infer the correct 
expression of the energy barrier from the accuracy of the 
data fit. On the other hand, the driving fields considered 
(H = 0.3, ...,0.6) are of the same order as the critical 
value (H c « 1.37) while Eq. (g|) is believed to be valid 
only in the limit H <C H c . 



V. DISCUSSION AND CONCLUSION 

We investigated numerically the creep motion of a 
driven interface of a RFIM model in the limit of low tem- 
peratures and small driving fields. We found that the in- 
terface velocity obeys an Arrhcnius law, which was inves- 
tigated in detail. We assumed that the prefactor of the 
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FIG. 9. On rescaling the interface velocities with the nu- 
merically determined energy barrier E(H) the data shown 
in Fig. M coalesce onto one single curve. As in the two di- 
mensional case this behavior shows that c(H) « const [see 
Eq. §]. 
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FIG. 10. Equation (^[) results in a logarithmic dependence 
of the energy barrier on the driving field for /i — > 0. The figure 
shows the interface velocities, which are rescaled according to 
Eq. §. 



Arrhenius law can be written as C(H,T) ~ c(H)T~ x . 
Applying this Ansatz to the numerically determined in- 
terface velocities, we find a positive exponent x for the 
two and three dimensional model. Additionally, our re- 
sults suggest that c{H) is independent of the driving 
field in both cases. These results are in contradiction 
to a renormalization group analysis [p"l|| in which (i) x is 
claimed to be negative and (ii) c(H) is found to exhibit 
a significant field dependence [c(H) ~ H 17 with a > 0]. 
In particular the opposite sign of the exponent x is re- 
markable. 

Knowing the prefactor C(H,T), it is possible to in- 
vestigate the driving field dependence of the energy bar- 
rier E(H). Our numerical results are in agreement with 
the assumption that the energy barrier depends logarith- 
mically in both dimensions on the driving field. Again 
this result is in contradiction with both phenomenologi- 
cal theories and renormalization group approaches, which 
conjecture an algebraic behavior |]TT|— |l3[| . The logarith- 
mic behavior can be explained if one assumes that the 
exponent of the algebraic behavior tends to zero. But 
in ]Tl|-p^] a finite value of the corresponding exponent is 
predicted. 

Thus our analyses reveal that the driven interface of a 
RFIM displays creep motion in the limit of low tempera- 
tures and small driving fields characterized by an Arrhe- 
nius law as predicted by phenomenological and renormal- 
ization group approaches. The details of the Arrhenius 
law (prefactor and energy barrier) differ, however, from 
the predicted behavior. Further investigations are needed 
to understand these differences. 
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